Bond Algebras and Exact Solvability of Hamiltonians: 
Spin £=1/2 Multilayer Systems and Other Curiosities 
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We introduce an algebraic methodology for designing exactly-solvable Lie model Hamiltonians. 
The idea consists in looking at the algebra generated by bond operators. We illustrate how this 
method can be applied to solve numerous problems of current interest in the context of topological 
quantum order. These include Kitaev's toric code and honeycomb models, a vector exchange model, 
and a Clifford 7 model on a triangular lattice. 
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I. INTRODUCTION 

Whenever one is interested in studying a new phys- 
ical phenomena whose effective model includes degrees 
of freedom (spins, fermions, bosons, etc.) which are 
strongly coupled, one attempts to invoke simplifying 
assumptions hoping that the resulting problem repre- 
sents the relevant minimal model for the phenomenon 
at hand. Finding exactly-solvable models is always wel- 
come. This paper is about a general methodology to gen- 
erate exactly-solvable Hamiltonians by concentrating on 
the algebra generated by algebraic objects called bonds. 
We have already used this methodology in Ref. [l[ , where 
we solved a doped orbital compass model in two and three 
space dimensions, although we did not explained the gen- 
erality of the mathematical approach. A goal of this pa- 
per is to present this methodology in full detail and show 
that by using the algebra of bond operators one can eas- 
ily construct whole families of exactly-solvable models, 
several of these displaying topological quantum order Q , 
and spin liquid behavior. 

For the sake of clarity, we will focus on quantum lat- 
tice systems which have N s = Jl^Li sites, with 
the number of sites along each spatial direction (i, and 
D the dimensionality of the lattice. The connectivity of 
the lattice and its general topology are of paramount im- 
portance. Associated with each lattice site i £ Z Ns there 
is a Hilbert space Hi of finite dimension T>i. The to- 
tal Hilbert space is the tensor product of the local state 
spaces, H = Hi, in the case of distinguishable subsys- 
tems (or a proper subspace in the case of indistinguish- 
able ones), and its dimension is T> = Yii=i ^i- 

Let us first start with an intuitive introduction to the 
key concept of bond algebras. We consider situations in 
which the Hamiltonian of a system H , whose state space 
is H, can be written as a sum of quasi-local terms or bonds 
{{h R }), 



H = ^ a R h R , 

R 



(1) 



general, the operators Hr will generate a certain (bond) 
algebra Q whose dimension is 0(T>). To simplify the de- 
scription, in the following, we are going to concentrate on 
semisimple Lie algebras. Notice that we do not constrain 
ourselves to a particular representation of the algebra. 

It may happen that the Hamiltonian itself is an element 
of a subalgcbra of Q of dimension polylogP. If such is 
the case, H, which represents a Generalized Mean-Field 
Hamiltonian (GMFH), is exactly-solvable Q, and there is 
a polynomially in log T> efficient algorithm to diagonalizc 
it 3. We say that H is exactly-solvable when an arbi- 
trarily chosen eigenvalue, and an appropriate description 
of the corresponding cigenstate, can be obtained and rep- 
resented to precision e by means of a classical algorithm 
efficient in log I? and 1/e. This definition, motivated by 
complexity theory, yields a sufficient criterion for exact- 
solvability. A particular case of exact solvability is when 
the spectrum can be expressed in closed form. 

A main contribution of this paper is to propose a 
methodology to generate such Hamiltonians by using two 
mathematical principles that will become evident in the 
following sections. In all cases, this methodology rests on 
(1) Topological constraints that are related to the con- 
nectivity of the lattice Hamiltonian or graph. In sev- 
eral instances, it further relies on the existence of (2) 
gauge symmetries. These symmetries allow a decompo- 
sition of the Hilbert space into sectors. The operators 
{h R } belong to the lowest dimensional representation of 
the algebra on the Hilbert space, or its sub-spaces. In 
the following we illustrate the bond algebra methodol- 
ogy by showing some tutorial examples of known trivially 
exactly-solvable problems. 



A. Ising model 

A simple example is afforded by the Ising model on a 
hypercubic lattice of N s sites, 



where an is a c-number (H must be an Hermitian oper- 
ator) and R includes a finite number of lattice sites i. In 



Rising — — ^ J<Ji<Jj. 



(2) 
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The bonds bij = OiOj satisfy a simple Ising (Abclian) 
type algebra defined on a T> = 2 Ns -dimensional space 
(the span of the original Ising system) : 

[6y-Az]=0, 6?. = 1, (3) 

since Oi = ±1. All classical Hamiltonians are extreme 
cases of GMFHs: Its spectra are trivially determined. 



B. Transverse field Ising chain 

The Hamiltonian of a single transverse field Ising chain 
of length N s reads 



-Htfim = -^{JiOiOi +1 + haf), 



(4) 



where erf = x, y, z) represent Pauli matrices. To make 
clear the algebraic connection that will follow, let us de- 
note the two terms (transverse field and bond variables) 
as follows: 



A? 



V V 

a l a r 



(5) 



In terms of these, the Hamiltonian of Eq.((4]) obviously 
reads 



H 



TFIM 



E^ 



(6) 



with interaction terms satisfying 



[Af 

[%,*j,k] 
{AD 2 







- 1 . i 



(7) 



= {A*,A i 

, i ^ j, k 

1 = fe) 2 , 

which forms an so(2jV s ) (polylogP) algebra with T> = 
2 N °. 

Note that the bond algebra encapsulated in the rela- 
tions above is invariant under the flip of any transverse 
field locally. The transformation 

Af - -At (8) 

effects hi — ► —hi at the lattice site i. Indeed, all that 
a flip of local fields does is to leave the spectrum un- 
altered while permuting the eigenstates amongst them- 
selves. In more conventional terms, the invariance of the 
spectrum mandated by the invariance of the bond alge- 
bra under the transformation of Eq. © is seen by noting 
that a similarity transformation with the local unitary 
(and Hermitian) operator [7* 



of sets 



= -of (9) 

while leaving af and thus Ai^ + i invariant. The spectrum 
of Eq.© can be determined by performing a Jordan- 
Wigner transformation to free fermions. Equivalently, 
it may noted that the bond algebra of a tight-binding 
spinlcss Fermi model (with pairing terms) is equivalent 
to that of Eq.©. 



C. Orbital compass chain model 

This model was introduced in Q . It consists of a D = 1 
dimensional system with alternating xx and yy interac- 
tions. Namely consider a chain of length N s in which 
the Hamiltonian is given by 



Hqc\ 



E 

i=l,3 5 5, 



J x 



E 

i=2,4,6,- 



Jy.iofcr. 



i+1 



Let us define the even and odd bonds by 



A, 



y y 

T 2m a 2m+1 



B„ 



'2m-l u 2m- 



(10) 



(11) 



They satisfy the following algebra (V = 2 s ) 

[A m ,A n ] = = [B. m , B n ] 

{A m ,B m } = Q = {A m , B m+ i}, 

[A m ,B n ] = , \m-n\ > 1 

{A m f = l = (B m ) 2 . 



(12) 



This algebra is identical to the algebra of bonds of 
Eqs.([7]). For an open chain, there are no boundary con- 
ditions on the bonds in either problem. If we enabled in- 
teractions (both exchange and transverse fields) on only 
one half of the chain (that is, if the sum in Eq.(d]) would 
extend, for even N s , only from 1 < i < N s /2) and add 
N s /2 non-interacting spins, then the number of inter- 
action terms in Eq.® and Eq. (fTT)|) . their algebras (and 
dimension of their representations), and the size of the 
Hilbcrt space on which both systems are defined are iden- 
tical. In that case, the partition functions are identi- 
cal up to a trivial multiplicative factor (after identifying 
J% and hi Jx.i) 



Z 0C u{N s ) = 2 N ^ 2 Z TFm (N s /2). 



(13) 



Such a relation was indeed found by [J] by an explicit di- 
agonalization of the Fermi bilinear found after a Jordan- 
Wigncr transformation performed on -ffoCM- Here we 
arrived at the same result by a trivial application of the 
methodology of bond algebras. 



D. Kitaev's Toric Code Model 

Kitaev's toric code model Q is defined on a square 
lattice with L x L = N s sites, where on each bond (or 
link) (ij) is an S = 1/2 degree of freedom indicated by a 



Pauli matrix er^ . The Hamiltonian acting on a T> = 2 
dimensional Hilbert space is 



2N S 



H 



A" 



E^-E 5 , 



(14) 



with Hermitian operators (whose eigenvalues are ±1) 



A„ = 



n 



Bp = 



n 



(15) 



(ij)<Estar(s) 



(y)Splaquette(p) 
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where B p and A s describe the plaquette (or face) and 
star (or vertex) operators associated with each plaquette 
p, and each site s of the square lattice. The reader may 
want to consult Refs. 0, [als! f° r notation purposes. 

That the D = 2 Kitaev's toric code model is identical 
to two decoupled Ising chains 0, 0] is immediately seen 
by looking at the bond algebra. The algebra of the bonds 
given by Eq. (fT"5|) is trivial, it is an Ising (Abelian) type 
algebra 

[A S ,B P ] = [A S ,A S ,] = [B p ,B p ,}=0, 
{A s f = 1 = [B p )\ (16) 

For periodic boundary conditions one has the additional 
constraint 



(17) 



It is very easy to realize that the Hamiltonian for two 
decoupled Ising chains, each of length N s 



Hi = - °~sO s +i 

8=1 



^2 t p t p+i 
i 



(18) 



with o~ s = ±1 and r p = ±1, displays an identical bond 
algebra to Eqs . (fT6|) . with the same representation. Thus, 
one can immediately write down the partition function 
0I| 



Z K = (2cosh/3) 2iv "(l + tanh Ws (3) 



(19) 



where (3 = l/(k,BT), and T is temperature. Moreover, 
the bond algebra of Kitaev's toric code model is identi- 
cal to that of Wen's plaquette model JTo] which proves 
the equivalence of the two systems [8j. It is worth- 
while to note that Eq. (p~9|) is also the outcome of a high- 
temperature series expansion (llj |. 

Thus, Kitaev's toric code model is identical to a one- 
dimensional Ising system. This statement has ramifica- 
tions for the stability of quantum memories — an item 
that we investigated in detail early on 0, This 
mapping allows not only an evaluation of the partition 
function but also a direct computation of all correlators in 
Kitaev's toric code model. For a detailed explanation see 
Refs. 0, H). In particular, see subsections (XIII A,B) as 
well as footnotes [61-63] of Ref. 0]. The equations of mo- 
tion with uncorrelatcd noise are insensitive to a change 
of basis. Consequently, the dynamics and thermal effects 
present in one-dimensional systems rear their head also 
in Kitaev's toric code model. In particular, the system 
is unstable to thermal noise — a phenomenon that we 



coined thermal fragility 




and has been recently 
151 1. Our bond algebraic 



uity 17 

confirmed by others [12|, LL3J, [1 
mapping enables an immediate extraction of crossovers 
for finite size systems fl6j |. 



E. Plaquette model in a transverse magnetic field 



This model [19( is defined on a square lattice as follows: 



Hr 



(20) 



where F, = o y l+&i o* + ~ x + ~ y o y l+&y , with e M representing 
unit vectors along the [i direction in the lattice. 

Setting, d = af, and Fi = F^ with i* = i + \{e x + 
e y ), and a diagonal chain coordinate j along the (1,1) 
direction, that alternates between i and i*, the algebra 
of the interaction terms (bonds) in Eq. (|20p is 



{Fj, G J+ i} 
(Fi) 2 



0=[G jt G r ] 
= {F j ,G i - 1 }, 







= 1 = (G 4 ) 2 . 



(21) 



For a system with open boundary conditions, the algebra 
of this system is none other than that of a stack of de- 
coupled transverse field Ising chains (see subsection II Bl 
and Eq.((7|) in particular) all oriented diagonally along 
the (1,1) direction. Setting, in Eq.©, Aij = A^, we see 
that the algebra and the dimension of the Hilbcrt space 
in both problems are identical. Indeed, a more elaborate 
treatment finds that this system is none other than that 
of a transverse field Ising model [l9| precisely as we find 
by examining the bond algebra. The partition function 
is, therefore, exactly the same as that of a transverse field 
Ising model. 

In the next sections, we illustrate the power of our 
method by reviewing several, more challenging, known 
examples of exactly-solvable models whose solutions can 
be immediately achieved in this way, and then we turn 
to new models that we introduced and solved using these 
tools. We start by discussing Kitaev's honeycomb model 
2C | and show that no enlargement of the Hilbert space 
20 j nor a direct Jordan- Wigner mapping [2l| is necces- 
sary to solve this model in a very short and direct man- 
ner. Next, we will turn to new models. The first of these 
is the vector exchange model which forms a simple ex- 
tension of Kitaev's honeycomb model. We will later on 
show that all these models have in common a Clifford 
algebraic structure on-site and an Abelian structure off- 
site. This defines a simple class of GMFH of the so(2iV) 
type 0]. More general Lie algebraic structures can be 
also realized. 



II. KITAEV'S HONEYCOMB MODEL 

A. Spectrum from bond algebras 

Kitaev's honeycomb lattice model [13, [HJ is a member 
of a family of models whose Hamiltonians are elements of 



4 



the so(27V s ) algebra, where N s is the number of vertices 
of the honeycomb lattice, i.e., it is a GMFH. 

The model is defined by the following 5=1/2 Hamil- 
tonian (Fig. [1]) 



H 



a: — bonds y— bonds z— bonds 



(22) 
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FIG. 1: Kitaev's model on a honeycomb lattice and three 
types of bonds. On each vertex there is an S = 1/2 degree 
of freedom indicated by a Pauli matrix There are two 
different types of vertices. Thick-colored contours represent 
arbitrary paths drawn on the lattice, e.g., from site i to j. 



Let us define the bond operators 



H = x,y,z, 



(23) 



where defines a bond (nearest-neighbors). We can 
assign a uniform direction to all bonds by choosing for all 
nearest neighbor sites i and j, the bond defined by 
corresponds to j — i = e M (and not j — i = — e M ). From 
the commutation relations for SU(2) spins it is clear that 
(A ijk i = 5 lk + S jk + Su + Sji) 



Ah 



= {l- 5tiV ){l-{-l)^)A%Al u (24) 



and, moreover, it is clear that (A^) 2 = 1. For an arbi- 
trary set of bond operators A^, the Lie algebra Q gen- 
erated is 0{T>). However, the set of bond operators that 
appear in i/if h , {A^}h k , forms a Lie subalgebra of Q, 
because of the particular lattice topology: Bond oper- 
ators that share a vertex anticommute, otherwise they 
commute. This subalgebra is precisely so(2_/V s ), and the 
Hamiltonian being an element of that subalgebra is an 
GMFH. 



One could, in principle, stop here and diagonalize the 
problem in a Hilbert space of dimension T> = 2 Na . How- 
ever, there is a further simplification in this problem. 
The simplification is related to exploting the existence of 
gauge symmetries. Consider the anyon charge [2(| oper- 
ators 



^= n 



A% = ofoZ<4oio» s <T S 



where h a defines a particular hexagonal plaquette (see 
Fig. [T|). The {Ih a } operators have eigenvalues J/ ta = ±1, 
and they satisfy the following relations 



[h a ,I h J = 
[h a ,At 3 ] = 



1. 



(25) 
the 



which implies that [Ih a , Hx h ] = 0. In other words 
product of bonds taken around the hexagon h a in a uni- 
form orientation (either clockwise or counter-clockwise) 
is a gauge symmetry. The N s /2 operators {4J satisfy 
the global constraint 



IF* 



i. 



(26) 



This allows us to decompose the 2?-dimcnsional Hilbert 
space Ti into 2 Ns ^ 2 ~ 1 orthogonal Hilbert subspaces Ti^, 
each of dimension dim(7i, ) ) = 2 Ns ' 2+1 



2 N 3 /2-l 

r,= l 



(27) 



Each Hilbert subspace TL^ is characterized by a particular 
set of eigenvalues {4 o } and projector 



P, 



N s /2 

n 

a=l 



t + I h J h 



N s /2 



(28) 



The algebra satisfied by the projected bond operators 
A^j = P rj A^-P rj is so(2iV s ), but acts on a Hilbert (carrier) 

subspace of dimension 2 Ns ^ 2+1 

To determine the spectrum in each subspace we look 
for an oscillator realization of the algebra 



2'//, ,<•,'•, , /i = x, y, z, 
in terms of Majorana fermions c,, which satisfy 



{Ci,Cj} 



„t 



(29) 



(30) 



The smallest representation of N s Majorana fermion 
modes (N s even) is that in a 2 Ns ^ 2 dimensional Hilbert 
space. For reasons that will become clear later on, in 
Eq. (f2"9"|) we will set r)ij = 1 on all bonds parallel to the "x" 
or "y" directions and allow 77^ = ±1 on all vertical bonds 
(those parallel to the "z" direction in which j — i = ±e z ). 
There remains an additional degeneracy factor of two 



5 



(the Hilbcrt space dimension is of size 2 Nb ^ 2+1 while the 
representation of the bonds is on a Hilbert space of size 

It is straightforward to show that the bilinear combina- 
tions of Majorana fermions satisfy the same algebra and 
also constraints as the algebra of {Alj}H Kh ■ Notice that 
any connected open string product of bonds becomes a 
bilinear in Majorana fermions [2l[ (see Fig. [T]) 



c. . _ 4P1 4M2 . . . 4 _ 9„--f 



n 

(i,3) 



,(31) 



It turns out that all open strings having the same 
end points i\,iL and with alternating fi's (e.g. 
x, z, y, z, y, x, z, x) can be of only 4 types 



(±l,±i) 2c n c lL+1 , 



(32) 



and form a polynomial in the number of vertices (or 
bonds) Lie algebra. The correspondence between the 
anyon charge sector {Ih a } and the set {%} is 



II 



Ih 



(33) 



{ij)eh a 



The set {h a } spans all fundamental hexagonal plaque- 
ttes from which all closed loops T can be uniquely con- 
structed. The number of plaquettes {h a } is given by 
half the number of sites N s /2 as is the number of verti- 
cal bonds on which we may asign one of the two phases 
corresponding to r\ij — ±1. With the identification of 
Eq. ([33)) . both sides of Eq. (|29[) satisfy the same set of 
algebraic relations in each of the 2 x 2 Ns / 2 dimensional 
Hilbcrt subspaces. 

This mapping allows us to immediately write down 
the spectrum in each sector of {Ih a } and to reproduce 
the results of [2(| HH without the need for introducing 
two Majorana fermions per spin and then projecting out 
one (as in [20]) nor writing expliciting a Jordan- Wigncr 
transformation between fermionic and the spin variables 
(as was done in pl|). 

In a given sector r\ = {fjij}, we have the Majorana 
fermion representation of the Hamiltonian, 



H Kh , n - 2i) j rjijJ i . 

(ij) 



CiCj 



(34) 



where Jij = J x if i and j arc separated by an "1" type 
bond. Similarly. = J VtZ if i and j are linked by a 
"y" or "z" type bond. Within the ground state sector 
(Ih a — 1 f° r au plaquettes h a ), we may set {rjij = 1} and 
obtain the quasi-particle spectrum, k = (k x ,k y ) pol. [2l|. 



Ar = 2 J x smk x + 2 J v sin k v . 



(35) 



Our mapping allows for a closed-form solution only 
for a reduced set of sectors (such as the ground state 



sector). For the rest, we still can compute each eigenvalue 
and eigenvector with polynomial in N s complexity by 
using the Jacobi method Thus, it is not simple to 
compute the partition function of the model with the 
same complexity. One can write down a formal solution, 
as wc discuss below, but it is not a closed-form analytical 
solution in terms of simple functions. 



B. Partition function 

Although many results appear on the zero temperature 
behavior of Kitaev's honeycomb model, there are very 
few results at finite temperatures. An exception is Ref. 
[23| which provides a finite temperature metric analysis of 
Kitaev's honeycomb model. Related results are discussed 
in ||. 

The partition function includes contributions from all 
sectors and reads M 



Z v = Trexp[-f3H Kh!r] ], 



(36) 



N s 



is the number of hexagonal plaquettes. 



where 

In terms of the original spins of Eq. ([22 

' (PH Kh ) 2n 



n=0 



(2n)! 



(37) 



or, equivalently, in terms of the Majorana representation 
of Eqs.dMEl 



Z,=Tr£ 



(2n)! 



(38) 



The reason why in Eq. (|3T[) we keep only keep the 
even powers of Hk h is the following. By time reversal 
symmetry (due to the trace over erf and (—erf)), at any 
given site i we must have an even power of erf. Similarly, 
in the Majorana fermion representation [Eq. (|38[) ]. the 
trace of Ci is zero. For any odd power of Hx h , there is 
in any term resulting from the expansion of exp[— /3HK h ] 
at least one site for which we have an odd power of erf 
(or Ci). All of these terms vanish. 

Let us first consider the Majorana representation and 
focus on Z v . We will later on rederive these results within 
the original spin representation of Eq. (|2"2"|) . The local 
assignments {rjij} effectively relate Jij in a general sector 
to that in the sector {rjij = 1} by the transformation 



'J ■ 



(39) 



We claim that if a particular bond (Jij) appears as an 
odd power in a given term then it will give rise to a 
vanishing contribution when it is traced over. The proof 
of this assertion is trivial: 



E 







(40) 



=±1 
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for all odd p. 

The same conclusion follows within the original spin 
representation of Eq. ([22]) which as wc show below leads 
to Eq. ([57)) . Let us mark all the bonds J^ that would 
additionally appear to an odd power in the expansion 
of Eq. ([57)) . We claim that there are several possible 
topologies: 

(i) three odd bonds (odd powers of Jy) meet at a com- 
mon vertex. 

That is, we can have 

J^J^JT (41) 

with odd pi a (a = j, k, I) and with all of the bonds that 
touch j, k and I appearing to an even power in the ex- 
pansion of Eq. ([37)) . 

(ii) Closed or open contours of odd powered bonds ap- 
pear. 

In case (i), the spins at sites j, k and I appear to an 
odd power (the power is just the sum of the powers of the 
bonds that have one of these points at their end). In case 
(ii) , if the contour is open then the spins at the endpoints 
of the open contour must appear to an odd power and 
therefore leads to a term that vanishes upon taking the 
trace. If the contour is closed it leads to none other than 
the anyon charge within the contour C, 

ic = n ih « ( 42 ) 

h a &C 

with Ih a the product of bonds along a hexagonal loop. 
Using the relation a^a v = ie^ K a K for jx ^ v. we find that 
each spin i that lies on C (i € C) leads to a contribution 
erf. For any odd power p, [erf ] p has a vanishing trace. 

We can similarly, have both (i) and (ii). It is readily 
seen that all odd powered bonds lead to situations with 
cither odd powers of the spins at the endpoints and/or to 
closed contours which also lead to vanishing contribution. 

Returning to the sum of Eqs. ([3"6l l3"8)h we now have 

Z = 2 Nh - N " +1 'E[2 1 }, (43) 

where [i| 

Z x = exp[-^Tr[Mi]] | det(l + exp( / 37V 1 ))| 1 / 2 (44) 

is the partition function of the system in which all rjij = 1 
(the anyon free system). In Eq. l|4"3"]). the operator E 
projects out of Z\ only those terms that are even (hence 
the symbol E) in all exchange constants {Jij}. The ma- 
trices Ni and Mi depend on the constants {Jij} as de- 
tailed in Q. The same conclusion follows from the ex- 
pansion Eq. ([37)) . 

We note that for the particular set of exchange con- 
stants Jj,i+e„ = Jfi, (nearest neighbor Jij), Eq. (|4"3")) is 
indeed equal to 

Z = 2 Nh+1 W[Z{\ (45) 



where E z projects out of Z\ only the terms that have all 
of the powers of Ji t i+e z being even. In order to cover all 
of the topological sectors Ih a = ±1, in each hexagon h a 
it suffices to allow rjij = ±1 on all of the vertical bonds 
(parallel to the z direction), and rjij = 1 on all other 
bonds (parallel to the x or y directions) [2l| . 

It is worth emphasizing that different topological sec- 
tors {Ih a } lead to different Z v (and thus to different spec- 
tra as they indeed must). It is only after performing the 
trace in Eq. (|38[) that the common even powered terms 
are pulled out. These terms are the same in all {rjij} 
assigments. 

III. VECTOR EXCHANGE LATTICE MODEL 

A. Motivation 

Consider the Lagrangian density of fermions coupled 
to a vector gauge field A a with a = 0, 1, 2, 3. In a U(l) 
theory, the Lagrangian density describing the minimal 
coupling of fermions to the gauge field is given by 

£mi„ = Hn a d a - l a A a )4>, (46) 

where j a — (■y a y are the Dirac matrices. Within the 
U(l) theory, A a= o is the scalar potential and j4 a =i,2,3 
are the spatial components of the vector potential A. 
The minimal coupling term of Eq. (|46[) is augmented by a 
gauge-only term {\F ab F ab with F ab = d a A b - d b A a ). In 
the electroweak theory (SU(2) x U{1)), the A a in Eq.® 
is replaced by (A a — V a ) with the weak parity breaking 
field V a . 

Although the Lagrangian is quadratic in the fermion 
fields ■0, it is definitely not a simple quadratic form that 
can be exactly integrated out. This is due to the linear 
coupling to A a . The theory contains both free quadratic 
terms (e.g., those in ijj alone) and terms of the form 
■07 a A Q '0. These terms give rise to interactions such as 
the lowest-order exchange term shown in Fig.©. The 
coupling to the A a gauge field gives rise to the usual 
Coulomb interaction between fermions. 

The lowest-order interaction terms are those formed 
by two vertices as above. The slanted lines depict the 
fermions (ip) while the horizontal wavy line represents the 
photon propagator (D ab )- the propagator for the fields 
A a . Integrating out the gauge field gives rise to the usual 
Coulomb exchange (depicted in Fig.©) 

j(xh a iP(x)D ab (x, y)ip(yh b ip(y) (47) 

with 

D ab (x,y) = (A a (x)A b (y)) (48) 

the Coulomb propagator. The same formalism albeit 
with more indices applies to other vector gauges (e.g., 
the electroweak one). In the non-relativistic limit, the 
density-density interaction (the ip(x)~f°ip(x)ip(y)~f ip(y)) 
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FIG. 2: Vector exchange between two fermions. To lowest 
(quadratic) order in the gauge field, a bilinear between the 
Dirac matrices: ip(x)j a ip(x)D a t,(x, y)4>(y)'y b '4'(y) results. For 
a vector exchange (photon) propagator D a b(x,y) set to be 
J x , y on a lattice, the resulting system is precisely of the form 
of the vector exchange model of Eq. (|50|l . 



piece becomes important. That is, the Doo propagator 
becomes dominant for non-relativistic particles. 

The lattice gauge action for the fermions resulting from 
integrating out the vector gauges A a is not usually inves- 
tigated in lattice gauge theory calculations. It is correct 
but this is not the standard point of departure for lattice 
gauge calculations. What is typically done is to write 
terms of the form 



ip(x)j a 4>(x + e a ) exp[i^4 



x,x-\-e a 



(49) 



with A x , x +e a the line-integrated lattice gauge (3+1) vec- 
tor potential between nearest-neighbor sites. 

What we do in the following affords another way of 
investigating general minimally coupled actions. When 
integrating out the A a fields, we generate precisely in- 
teractions of the 77 type of Ea.(|I7]) with Eq.(|3g|). In 
what follows, we will investigate a simple lattice rendi- 
tion of such vector exchange system in which we set an 
exchange coupling between 7 matrices to be of amplitude 
Jij = D(i,j) with i and j denoting lattice sites. 



B. Exact solution of the vector exchange model 

A simple square lattice model that captures the 
fermionic vector exchange is given by 



H = J2 J ^a,iJa,j- 

(ij) 



(50) 



The geometry of the lattice is shown in Fig.©. The 7 
matrix index a for a given bond in Eq. (|50[) is fixed by 



the sites i and j. Here and throughout, we set the lattice 
constant to be one. 



IHf^ ^ ^ ^llllllllllllllllllllllll^ ^ ^ ^fttl 



ittf^^ ^iiiiiiiiiiiiiiiiiiiiiiii^ 



If"" 



t m 

72,i72,mfi 



t 1 1 



T4,i74j' 

-■ °j 



,t 



7i,i7i,jfc 



FIG. 3: The bonds in the system of Eq. (|50|l . At each ver- 
tex i there are four different types of bonds corresponding to 
7a,i7a,j interactions with different a — 1,2,3,4. (There are 
four different types of vertices.) The algebra of the bonds is 
trivial: bonds that share a site anticommute, disjoint bonds 
commute, and the square of any bond is one. Consequently 
the spectrum of the model can be immediately determined. 

The 7 matrices satisfy the algebra 

{la,i,lb,i\ = 2<5 q6 , [7a,i,7b,j] = 0, i ^ j. 



(51) 



The Hilbert space on which H acts on is, for a lattice 
of N s sites, of dimension 4 Na . The algebra of the bonds 
la,ila,j is familiar: it has the same simple characteristics 
of the bond algebra in Kitaev's honeycomb model. These 
algebraic relations do not change on projection to a state 
of fixed topological charge sector 



Pod )Pt 



(52) 
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We define the projector Pq to a topological sector by 

1 + I a Ia 



with 



Pa 



(53) 



(54) 



As in Eq. ([28"|) . /□ arc c- numbers: /□ = ±1, and prod- 
ucts of bonds around a plaquette Jp commute with 
the Hamiltonian ([Jp,if] = 0), and amongst themselves 
d/dj-fc'] = 0). The operators of Eqs. (f54")) constitute lo- 
cal (gauge) symmetries. The origin of the commutation 
relations is that at each vertex we have bonds of differ- 
ent 7 matrix flavors. Similar to the situation in Kitaev's 
honeycomb model, all bonds commute with the anyon 
charge operators of Eq. (j54")) , and 



(55) 



1^ = 1. 



The gauge symmetries {Ins} allow decomposition of the 
total Hilbcrt space into orthogonal subspaces of equal 
dimensionality. We divide the Hilbert space into equal 
sectors spanned by {/□}. There are 2 JVs_1 such sectors 
as the eigenvalues of Jp, for each of the N s plaquettes □, 
can attain one of two values (±1), and satisfy only one 
global constraint on a torus 



t. 



(56) 



As {ip} are good quantum numbers, we may diago- 
nalize the Hamiltonian in a Hilbert space of dimension 
4 JV S/ / 2 iV 5 -i = 2 N <+ 1 . Similar to our solution of the Ki- 
taev's honeycomb model, we may then work with the rep- 
resentation of the bonds as the product of two fermions. 

Within each anyon charge sector, the Hamiltonian is 
of the form of Eq. (f34"]) but on different size spaces. We 
can now introduce N s spinless fermion variables {d{\ on 
the Hilbcrt space of size 2^" by setting the bonds to be 
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v -iidi + d^idj + d}). 



(57) 



We thus arrive at a Fermi bilinear that is trivially diag- 
onalizable 



H = iY^VijJij(di + d\)(dj 

(ij) 



A). 



(58) 



In Eq. (|58p . we maintain the directionality that we em- 
ployed throughout in constructing the bond algebra in 
the case of Kitaev's honeycomb model: j — i = e x or 
e y . A trivial but important feature of Eo. (f5"8"]) is that the 
spectrum is symmetric about zero. This is so as there is 
a symmetry Aij — > — Aij in the representation chosen in 
Eq.([57]). 

The dimension of the Hilbert space is the same as that 
of the product of all plaquette charges ip = Il<tj)en Vij 



(there are (N s — 1) such Ising type operators with eigeval- 
ues ±1 leading to 2 Ns ~ 1 sectors multiplied by the size 
of the Hilbcrt space spanned by the N s fermions in a 
space of size 2^ multiplied by a degeneracy factor of 
two. Fixed anyon charges enable 2^= possible configura- 
tions (redundant) of r\ij that give rise to the same original 
Hamiltonian when projected onto a sector of fixed {ip}. 
These configurations of r\ij are related to each other by 
local Ising gauge transformations on the square lattice. 
That is, with arbitrary n = ±1 at any lattice site i, the 
local gauge transformation 



(59) 



leaves /□ invariant. 

Written longhand, the 4^= dimensional Hilbert space 
spanned by the 7 matrices decomposes as follows, 

4 JV a = p^- 1 ( num ber of sectors {J n }) 

x 2^" (Hilbert space spanned by N B fermions) 
x 2(rcmaining degeneracy of each state)]. (60) 

There is a degeneracy factor of 2 = A N -/{2 N -~ 1 x 2^) 
that remains after invoking the representation of Eq. (|57p 
in the space of size 2^. This is similar to the degener- 
acy factor of two in subsection III Al The Hamiltonian 
of Eq. (|58p is nothing but a tight- binding Hamiltonian 
agumented by pairing terms (an element of the so(2A^ s ) 
algebra) on which we may apply a Bogoliubov transfor- 
mation similar to [2l| which was defined on the square 
lattice. The solution to Eq. I|58p can be immediately writ- 
ten down. For equal to J x or J y for sites i and j sep- 
arated by one lattice constant along the x or y directions 
respectively, i.e. J io = {Jx5\i m -j x \,i + Jy^\i v -j s \,i), m tnc 
sector rjij = 1 (corresponding to the sector 7p = 1), wc 
have on Fourier transforming, 



~Hn d % d v 



(61) 



with qj: = \J x e * + J y e %k v] and = 2(J x sm.k x + 
J y sin ky). A Bogoliubov transformation gives the quasi- 
particle spectrum 

E n = t ±2 Pi , (62) 
with the zero eigenvalue being doubly degenerate. 

C. String correlation functions from symmetries 

The only correlators that can obtain a finite expec- 
tation value at zero and finite temperatures must, by 
Elitzur's theorem (26|, be invariant under all local sym- 
metries. In [2l|, this only allowed for string correlators 
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of the form of Eg. (|3~l]) to have non- vanishing expecta- 
tion values at finite and zero temperatures. These sym- 
metry only conditions did not need to invoke the Majo- 
rana fermion representation used in [27l ]. The consider- 
ations of [2l[ for Kitaev's model can be replicated mu- 
tatits muntandis for the vector exchange model. The 
local symmetries of Eq. (|54p . allow only for string corre- 
lators (whether open or closed) to attain a finite expec- 
tation value. This applies for both the ground state con- 
figuration as well as the more physically pertinent case 
of all finite temperatures, T > 0. The zero temperature 
selection rule was also noted by [25| . Closed loops of the 
form of Eq. (|3~lj) that span the entire system correspond 
to additional symmetries of the system. Relying on simi- 
lar symmetry conditions as in [fl , it is seen that the only 
orders that can exist are of a non-local nature. 



D. Isomorphic spin models 

We now discuss two spin representations of our exactly- 
solvable vector exchange model: 

(i) a spin S — 1/2 variant on a square lattice bilayer 
and 

(ii) a spin S = 3/2 on a square lattice. This latter 
variant was very recently also discussed by [24], [H| . 
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FIG. 4: A square lattice bilayer of a S — 1/2 system that 
represents Eg . (|50fl . Inserting Eq. (|64[) , we find the interactions 
schematically depicted above. The two layers are indexed by 
the first of the subscripts (1 or 2) in Eq. (|64[) . 



|25| for their J5 = 0. The spin liquid character of these 
systems is inherited from Kitaev's honeycomb model. 



1. A spin 1/2 model on a square lattice bilayer 

A possible representation of Eq. (fSTJjl is obtained at by 
setting 



7m 



Hcr M ) 
(tcr^ 1 ) 



74 



1 2 
-l a 



(63) 



with /j, = 1,2,3, cr M the Pauli matrices, and I2 the 2- 
dimensional unit matrix. The 7 matrices can arise from 
the tensor product of two S — 1/2 spins at each lattice 
site. 



S 3 
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2 
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2 





1 








1 





%/3 
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V 





2 


J 



(65) 



7M 


= °IA 




72, i 


V V 

= °l,i°2,i> 




73, i 


V z 

= ai,i a h> 






= °h- 


(64) 



A possible lattice topology realization for this is that of 
a square lattice bilayer shown in Fig. The subscript 
a in <7^ i is the bi-layer index. Inserting Eq. (|64[) into 
Eq. ([50|) leads to a spin S =1/2 Hamiltonian on a square 
lattice bilayer. 



gy 



and 



/ 

2 ' 



V 



,■ \/3 
1 2 


i 




\ 

-i 
-a^i 
&i / 



(66) 



2. A spin 3/2 model on a square lattice 

Another representation of the model of Eq. ([5D|) is in 
terms of S = 3/2 spins that reproduces the results of [24|, 



s'- 







\ 




\o -I J 



1% 

i 

0-: 



(67) 
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In terms of two spins of size S = 1 /2 
S* 



s z 

(S x ) 2 

(S y ) 2 
(S z ) 2 



1 z z 5 

— <7i (To H ■ 

2 1 2 4 



V3 



1"2 



°1 CT 2 



5 
4' 



{S y ,S z } = >/3oM. 
We can represent the 7 matrices by 



(68) 



7i 
72 
73 
74 
75 



r* = i=(S*) 2 -(S*0 2 ), 



1 



'71727374 = of erf 



(S 2 



5 
4' 



(69) 



IV. CLIFFORD ALGEBRAIC MODELS 



The commonality of all these exactly-solvable models is 
the presence of degrees of freedom that satisfy a Clifford 
algebra on-site 



{la.i,lb,i\ = 2<5 ab , 

and a commutative algebra off-site 

[la.i,lb,j] = , 1 ^ j, 



(70) 



(71) 



with a,b = ].,••• ,p. The exactly-solvable Hamiltoni- 
ans are then written as linear combinations of quadratic 
products of these 7 matrices. Regardless of the dimension 
of the representation of the 7 matrices, the Hamiltonian 
is always an element of so(2iV s )(in the examples worked 
out in this paper), and thus a GMFH. 

From the viewpoint of lattice connectivity, notice a 
fundamental difference between Kitaev's honeycomb and 
the 77 (vector-exchange) models. The coordination of 
the honeycomb lattice is Zh = 3, while the one for the 
77 lattice is z 7 = 4. This is the reason why one needs 
p = 3 anticommuting (Pauli) matrices in the first case, 



while p = 4 anticommuting (7) matrices are needed in the 
second model. The relation between p and the dimension 
of the matrix representation of 7 is the following: When 
p = 2q or p = 2q + 1, the matrix representation can 
be of dimension 2 q . This is the reason why Kitaev used 
Pauli matrices (q — 1) in his honeycomb model, while we 
had used Dirac matrices (q = 2) in the vector-exchange 
model. 

It is indeed obvious how to generalize these ideas to 
generate new exactly-solvable models of the so(2iV) type 
in arbitrary dimensions and for arbitrary lattice coor- 
dination. The idea consists in writing bond operators 
which are quadratic products of Clifford operators which 
are anticommuting on the same lattice site. The cardinal 
p of that set of operators will define the z of the lattice (its 
connectivity). For instance, suppose we want to have a 
lattice with z = 5. Then, a Shastry-Suthcrland-like con- 
nectivity lattice will do the job [28|]. Now, write down 
a Hamiltonian which is a linear combination of bilinears 
of 5 anticommuting 7 matrices that act upon a Hilbcrt 
space of dimension A Na . This model will be exactly solv- 
able. In this way, we can construct a new model in a 
cubic z = 6 lattice with p = 6 anticommuting matrices, 
or as we show now a triangular z — 6 lattice model with 
7 matrices of dimension 2 3 x 2 3 . 

Consider the p = 6 7 matrices 



7i, i 


= °IA 


72,i 


V V 

= °IA 


73,i 


= <,< 4 


74,; 


= a l,i a 3,i 


75,t 


z V 

= °\,i a h 


76,i 





(72) 



which form an on-site Clifford algebra. The model Hamil- 
tonian 



H 



(ij) 



7a,i7a,i 1 



(73) 



whose lattice geometry is shown in Fig.©, is exactly 
solvable. This model can also represent a tri-layer system 
with plaquette interactions. 

One can indeed realize that there is nothing special 
about the Lie algebra so(2iV). One can consider mod- 
els whose bond algebra forms any other semisimple Lie 
algebra, such as so(2iV + 1) where there are non-linear 
Bogoliubov transformations that diagonalize the prob- 
lem. It is important, though, that the number of gauge 
symmetries is enough to allow for a simple oscillator re- 
alization of the bonds. Otherwise, there is always the 
possibility to use the Jacobi algorithm 0] to numerically 
diagonalize the problem. 
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FIG. 5: A triangular lattice Clifford model that is exactly 
solvable. In each vertex there is a degree of freedom of di- 
mension 2 3 (and a Hilbert space of dimension 2 3JVs ). There 
are four different types of vertices and six different types of 
bonds 7„,i7oj. 



V. CONCLUSIONS 

The thesis of the current work is that even though 
the solution to many problems is hard, by disregarding 
the explicit microscopic degrees of freedom and focusing 
solely on the algebraic relations that the bond variables 
satisfy in a Hilbert space of a fixed dimensionality wc 
may map one initially seemingly hard problem onto an- 
other problem whose solution is easier. This mapping 
does not rely on explicit real space forms for the trans- 
formations (although these can be written down in some 
cases). Nor does it rely on enlarging the Hilbert space 
and then making a projection onto a physical sector as in 
[20| | . It is important to emphasize, though, that the di- 
mension of the representation of the algebra is of crucial 
importance - not only the algebra and set of constraints 
itself. 

The explicit real space mappings - no matter how com- 



plicated their forms are - are irrelevant. The spectra and 
all non- vanishing correlators may be determined from the 
algebra alone. [See @ for a derivation of all correlators 
in Kitaev's Toric code model by this method.] In the 
current work, we illustrated how the energy spectra may 
be determined. The partition function and the density 
of states associated with the spectrum arc related by a 
Laplace transform. It may be easily seen also from the 
partition functions themselves that if two systems display 
the same bond algebra on a space of the same represen- 
tation then their spectra are identical. 

In the case of generic GMFHs the Jacobi method al- 
ways enables a solution of its spectrum with polynomial 
complexity ||. There are situations, like the model ex- 
amples presented in this paper, where we can determine 
the spectrum of certain sectors (Hilbert subspaces) in 
closed form, i.e., by quadrature. This will happen when- 
ever the effective matrices that need to be diagonalizcd 
have dimension smaller or equal to 4 x4. The decom- 
position of the Hilbert space into these individual decou- 
pled subspaces is rooted in the existence of local ( gauge ) 
symmetries. Within each subspace, there is an oscilla- 
tor realization (e.g., in terms of Majorana fermions) of 
the bond algebra. It is important to emphasize that ex- 
act solvability does not imply that we can compute the 
density of states, and thus the partition function, with 
polynomial complexity Kitaev's honeycomb model is an 
example of a system whose energy eigenvalues can be 
determined with polynomial complexity but whose total 
partition function cannot since its density of states is not 
determined with the same complexity. By contrast, Ki- 
taev's Toric code model constitutes an example of a sys- 
tem where not only the spectrum but also its partition 
function is exactly solvable 0, HI ■ 

Note added in proof. This work (and our vector- 
exchange model) was conceived in 2007. A physical 
model whose exact solution was enabled by our bond 
algebra mapping is detailed in [l[. During the time in 
which the current work was summarized, three works ap- 
peared [H, [2f| [29[ that introduce and solve variants and 
exact forms of the vector-exchange model that we intro- 
duced here. In particular, (2{| presented a S = 3/2 spin 
system and the vector-exchange model on the Shastry- 
Sutherland [28| decorated square lattice. 
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